Universality in metallic nanocohesion: a quantum chaos approach 
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Convergent semiclassical trace formulae for the density of states and cohesive force of a narrow 
constriction in an electron gas, whose classical motion is either chaotic or integrable, are derived. It 
is shown that mode quantization in a metallic point contact or nanowire leads to universal oscillations 
in its cohesive force: the amplitude of the oscillations depends only on a dimensionless quantum 
parameter describing the crossover from chaotic to integrable motion, and is of order 1 nano-Newton, 
in agreement with recent experiments. Interestingly, quantum tunneling is shown to be described 
quantitatively in terms of the instability of the classical periodic orbits. 



An intriguing question posed by Kac [1] is, "Can one 
hear the shape of a drum?" That is, given the spec- 
trum of the wave equation [1] or Schrodinger's equa- 
tion for free particles [2] on a domain, can one infer 
the domain's shape? This question was answered in the 
negative [1,2]; nevertheless there is an intimate relation 
between the two. In the context of metallic nanoco- 
hesion [3-10], a related question has recently emerged: 
"Can one feel the shape of a metallic nanocontact?" It 
was shown experimentally [3] that the cohesive force of 
Au nanocontacts exhibits mesoscopic oscillations on the 
nano-Newton scale, which are synchronized with steps of 
order 2e 2 jh in the contact conductance. In a previous 
article [4], it was argued that these mesoscopic force os- 
cillations, like the corresponding conductance steps [11], 
can be understood by considering the nanocontact as a 
waveguide for the conduction electrons (which are re- 
sponsible for both conduction and cohesion in simple 
metals). Each quantized mode transmitted through the 
contact contributes 2e 2 /h to the conductance [11] and a 
force of order ef/Xf to the cohesion, where Xf is the de 
Broglie wavelength at the Fermi energy sf- It was shown 
by comparing various geometries [4] that the force oscil- 
lations were determined by the area and symmetry of 
the narrowest cross-section of the contact, and depended 
only weakly on other aspects of the geometry. Subse- 
quent studies confirmed this observation, both for generic 
geometries [5,7,8,10], whose classical dynamics is chaotic, 
and for special geometries [6,9], whose classical dynamics 
is integrable. The insensitivity of the force oscillations to 
the details of the geometry, along with the approximate 
independence of their r.m.s. size on the contact area, was 
termed universality in Ref. [4]. A fundamental explana- 
tion of the universality observed in both the model cal- 
culations [4-10] and the experiments [3] has so far been 
lacking. 

In this Letter, we derive semiclassical trace formu- 
lae for the force and charge oscillations of a metallic 



nanocontact, modeled as a constriction in an electron 
gas with hard- wall boundary conditions (see Fig. 1 inset), 
by adapting methods from quantum chaos [12-16] to de- 
scribe the quantum mechanics of such an open system. 
It is found that Gutzwiller-type trace formulae [12-16], 
which typically do not converge for closed systems, not 
only converge, but give quantitatively accurate results for 
open quantum mechanical systems, which are typically 
more difficult to treat than closed systems by other meth- 
ods. Using these techniques, we demonstrate analytically 
that the force oscillations SF of a narrow constriction in 
a three-dimensional (3D) electron gas (i) depend only on 
the diameter D* and radius of curvature R of the neck, 
(ii) have an r.m.s. value which is independent of the con- 
ductance G of the contact and depends only on a scaling 
parameter a which describes the crossover from chaotic 
to integrable motion, and (iii) are proportional to the 
charge oscillations induced on the contact by the quan- 
tum confinement. Furthermore, we show (iv) that quan- 
tum tunneling through the constriction is determined by 
the instability of the classical periodic orbits within the 
constriction, and that the force and charge oscillations 
are suppressed only weakly (algebraically) by tunneling, 
unlike conductance quantization, which is suppressed ex- 
ponentially [11]. Conclusion (ii) is specific to 3D con- 
tacts, and breaks down for, e.g., two-dimensional (2D) 
nanowires, where rms SF cx G -1 / 2 . Conclusions (i), (ii), 
and (iv) are unchanged when electron-electron interac- 
tions are included within the Hartree approximation. 

The properties of simple metals are determined largely 
by the conduction electrons, the simplest model of which 
is a free-electron gas confined within the surface of the 
metal. Here we take the confinement potential to be a 
hard wall; the effects of interest to us are virtually un- 
changed when one considers a more realistic confinement 
potential [7]. The grand canonical potential £1 is the 
appropriate thermodynamic potential describing the en- 
ergetics of the electron gas in the nanocontact [4], and 
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where g(E) is the electronic density of states (DOS) and 
(3 is the inverse temperature [17]. The total number of 
electrons in the system is 



N-= J dEf(E)g(E), 



(2) 



where f(E) is the Fermi-Dirac distribution function. The 
DOS of an open quantum system, such as that shown in 
Fig. 1 (inset), is given in terms of the electronic scattering 
matrix S(E) by [18] g(E) = (2iri)- 1 Tr{S^(E)dS/dE - 
H.c.}, where a factor of 2 for spin has been included. 

The DOS can be decomposed [13,14] in terms of a 
smooth Weyl contribution g(E) and a fluctuating term 
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where k E — (2mE/h ) x l 2 , V is the volume of the electron 
gas, S is its surface area, and C = \ J da (l/i?i + I/-R2) 
is the mean curvature of its surface, R1.2 being the prin- 
cipal radii of curvature. The first three terms in Eq. 
(3) are macroscopic, while Sg determines the mesoscopic 
fluctuations of the equilibrium properties of the system. 
Inserting Eq. (3) into Eqs. (1) and (2), and taking the 
limit of zero temperature, one finds 
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where kp = 2tt/\f is the Fermi wavevector. The correc- 
tions to Eqs. (4) and (5) at finite temperature may be 
evaluated straightforwardly [13], and are quite small at 
room temperature, since Sp/k E > 10 4 K. 

The cohesive force of the nanocontact is given by the 
derivative of the grand canonical potential with respect 
to the elongation, F = —dft/dL. Under elongation, the 
contact narrows and its surface area S increases. The 
increase of S under elongation would lead to a macro- 
scopic surface charge by Eq. (5). This is due to the 
hard-wall boundary condition, which leads to a deple- 
tion of negative charge in a layer of thickness <~ Xp at 
the boundary [19]. The macroscopic incompressibility of 
the electron gas can be included by imposing the con- 
straint iV_ = const. [20], where N_ is given by the first 
three terms in Eq. (5). The macroscopic electronic charge 
— eiV_ is neutralized by the equal and opposite positive 
charge of the jellium background. The net charge imbal- 
ance on the nanocontact (neglecting screening) is thus 



SQq = —eSN-, which we will show to be quite small — on 
the order of a single electron charge. Differentiating Eq. 
(4) with respect to L with the constraint iV_ = const., 
one finds 



F = - — 
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where cry = £Fkp /16n is the surface energy of a noninter- 
acting electron gas [4] at fixed V and AF top = 46^/9 Af- 
The reduction of the surface energy by a factor of 5 has 
been discussed by Lang [19]. The second term on the 
right-hand-side of Eq. (6), termed the "topological force" 
by Hopplcr and Zwerger [5] since it depends only on the 
topology of the cross-section in the adiabatic limit, is re- 
duced by a factor of 2.5. Importantly, since the constraint 
7V_ = const, differs from the constraint V = const, 
used in previous work [4-6,8-10] only by terms of or- 
der (kpD*)^ 1 , the mesoscopic fluctuations SF and 5N- 
are quite insensitive to the choice of constraint. 

The fluctuating part of the DOS 5g may be evaluated 
in the semiclassical (stationary-phase) approximation as 
a sum over the periodic classical orbits of the system 
[12-16]. For closed systems, the sum over periodic orbits 
is generically not convergent, and a broadening of the en- 
ergy structure in 5g(E) must be introduced by hand [13]. 
However, we shall see that for an open system, such as a 
nanocontact, the periodic orbit sum converges; the finite 
dwell-time of a particle in an open system introduces a 
natural energy broadening. 

Let us first consider the case of a 2D nanocontact. For 
a finite radius of curvature R, there is only one unstable 
periodic classical orbit (plus harmonics), which moves up 
and down at the narrowest point of the neck. One obtains 
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where the Lyapunov exponent \ 01 the primitive periodic 
orbit satisfies exp(x) = 1 + D*/R + + D*/R) 2 - 1. 
Eq. (7) diverges when x — *• 0, i.e., when R — > 00. In that 
limit, the nanocontact acquires translational symmetry 
along the z axis, so that a generalization of the Gutzwillcr 
formula obtained by Creagh and Littlejohn [15] must be 
used, which gives a finite result. In this limit, the motion 
is classically integrable. One can treat small deviations 
from translational symmetry via perturbation theory in 
1/R. The resulting asymptotic behavior for large R may 
be combined with the result [Eq. (7)] valid for small R 
to construct the following interpolation formula, valid for 
arbitrary R: 
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where C(x,y) = cos(x)C(y) — sin(x)S(y), with C and S 
Fresnel integrals. In Eq. (8), the specific shape of the 
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nanocontact was taken to be D(z) = D* + z 2 /R. For 
a discussion of related interpolation formulae, see Ref. 
[16]. Classically, only the case R = oo is integrable. 
But semiclassically, there is a smooth crossover between 
the strongly chaotic (R — > 0) and the nearly integrable 
(R — > oo) regimes. The scaling parameter describing this 
crossover is 



a = Lj \J\pR. 



(9) 



We refer to a as the quantum chaos parameter, since the 
quantum fluctuations of the system correspond to those 
of a chaotic system when a ^> 1 and correspond to those 
of a quasi-intcgrable system when a<l. 

Fig. 1 shows a comparison of the semiclassical result 
flsc = 9 + &7hrt an d a numerical calculation of g using 
a recursive Green's function technique [10]. The agree- 
ment of the semiclassical result and the numerical calcu- 
lation is quite good, even in the extreme quantum limit 
G < 2e 2 /h. The small discrepancy is of the size expected 
due to diffractive corrections [13] from the sharp corners 
present in the geometry studied numerically, where the 
nanocontact was connected to straight wires of width 
kpD = 52 for technical reasons. 

The denominator sinhnx in Eqs. (7) and (8) describes 
the effects of tunneling. In the limit R ^> D* , the Lya- 
punov exponent \ ~ * \J2D* / R, and one recovers the 
WKB approximation of Ref. [4]. In the opposite limit 
R <C D* , sinhx — ► D*/R, so Sg is suppressed relative 
to the value expected in the WKB approximation (which 
neglects tunneling) by a factor of \/2R/D*. In the adia- 




FIG. 1. Inset: Schematic diagram of a metallic nanocon- 
tact. (a) Conductance G and (b) DOS <?(£f) for 2D nanocon- 
tacts with a m 5 versus the contact diameter D* . g is normal- 
ized to the area A of the region. Solid curves: semiclassical 
result based on Eq. (8); crosses with error bars: numerical 
results using the method of Ref. [10]. Lower curves in (a) and 
(b): R « \f; upper curves (offset vertically): R m 170Af- 



batic approximation, the energies of the transverse modes 
in the point contact are e n (z) — (h 2 /2m)(irn/ D(z)) 2 = 
e„(0) — muj 2 z 2 /2 + ■■■ and the probability that an 
electron of energy E in mode n will be trans- 
mitted through the point contact is [11] T n (E) ~ 
(1 + exp{-27r[£; - e n (0)]/7ici; n }) _1 . The quality of the 
conductance quantization thus decreases exponentially 
with the parameter ?kj n /As n ~ ir^ 1 y/2D* / R, where 
Ae„ = e n — e„_i, while the DOS fluctuations Sg are 
suppressed only inversely proportional to this parameter. 
The fact that the suppression in each case depends only 
on the ratio D* / R implies that the suppression of Sg, like 
the degradation of conductance quantization, is a conse- 
quence of tunneling. Indeed, it is the rounding of the 
DOS due to tunneling that causes the sums over n in 
Eqs. (7) and (8) to converge. That quantum tunneling 
through a point contact can be expressed purely in terms 
of the instability \ 01 tne classical periodic orbits within 
the contact is remarkable. 

Let us now consider the experimentally relevant case 
of an axially-symmetric 3D nanocontact. For finite R, all 
classical periodic orbits lie in the plane of the narrowest 
cross section of the contact; however there are now count- 
ably many distinct families of singly-degenerate periodic 
orbits [13,9], labeled by their winding number w about 
the axis of symmetry z and by the number of vertices 
v > 2w. The interpolation formula for Sg, describing the 
crossover from the chaotic regime a > 1 to the integrable 
regime a <C 1, is 
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xC(k E L vw - ^,a^vs\n(t> vw k E /k F ), (10) 
where <j> vw — nw/v, f vw = 1 + 9(v — 2w), and 



exp(xvw) = 1 + 



with L vw — vD* sin <p vw the length of a periodic orbit. 
We emphasize that the double sum over w and v in Eq. 
(10) converges due to the finite Lyapunov exponent Xvw 
In Eq. (10), higher-order terms in the small parameter 
1/fcp-D* ( < 0.21 for contacts of nonzero conductance) 
have been omitted. 

The mesoscopic force and charge fluctuations are cal- 
culated by inserting Eq. (10) into Eqs. (1), (2) and (6). 
In order to demonstrate the universality of the force os- 
cillations, it is necessary to make some physically rea- 
sonable assumptions regarding the scaling of the geome- 
try when the nanowire is elongated. It is natural to as- 
sume that the deformation occurs predominantly in the 
narrowest section, where the wire is weakest. This as- 
sumption, combined with the constraint of incompress- 
ibility iV_ = const., implies D* 2 L w const. Furthermore, 
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the radius of curvature R oc L 2 /(D — D*), where D 
is the diameter at ±L/2, which implies dhiR/ dh\L = 
2 + (dlnD*/dlnL)/(D/D* - 1) w 2. Thus the quantum 
chaos parameter a « const, under elongation. 

Using these assumptions about the scaling of the ge- 
ometry with elongation, the derivative with respect to L 
in Eq. (6) can be evaluated; the general formula for SF 
is rather lengthy, and will be presented elsewhere. Here 
we give only the limiting behavior of the leading-order 
semiclassical results: 



SF 
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a^>l w—1 v—2w 



f vw sin(k F L 



JF.-^VV; %sin(fc F ^-3W2), (12) 
Af ^-"i •^ _ r / f 

a<g;i u>— 1 f— 2uj 

where 6„ = 3vn/2 — 7r/4. is an oscillatory function of 
kpD*; the conductance of the contact is also determined 
by k F D*, indicating that the force oscillations are syn- 
chronized with the conductance steps, as shown in Ref. 
[4] and observed experimentally [3] . 

The rms amplitude of the force oscillations may be 
readily calculated from Eqs. (11) and (12). We find that 
rms 5F is independent of D* , and, apart from small cor- 
rections due to tunneling when R <C D* , depends only 
on the quantum chaos parameter a: 



rms SF = 



0.36208a" 1 a > 1, 



0.58621-^, a<l. 



(13) 



The result for a«l agrees with the result for a straight 
wire (a = 0) derived previously by floppier and Zwerger 
[9]. Eq. (13) is also consistent with previous results based 
on the WKB approximation [4]. For a realistic geome- 
try of the nanowire [3], one expects both the radius of 
curvature and the elongation to be on the scale of Af, 
implying a ~ 1. There is also experimental evidence [21] 
of exceptional geometries with R 3> Af, implying a<Cl. 
Thus the mesoscopic oscillations of the cohesive force are 
expected to be universal rms SF ~ Bf/Xf ~ InN in 
monovalent metals, in agreement with all available ex- 
perimental data [3]. 

In nanowires lacking axial symmetry, e.g., with an as- 
pect ratio a » 1, one can show that rms SF <~ ae F /\ F . 
However, such shapes are energetically highly unfavor- 
able due to the increased surface energy. Eq. (13) is 
therefore expected to describe all spontaneously occur- 
ing nanocontacts. 

Eq. (10) and the assumption D* 2 L = const, imply 
that the force and charge oscillations are proportional 
to each other in 3D nanocontacts: SF = —e F SN-/L + 
0{l/k F D*). In an interacting system, the charge os- 
cillations are screened [8], and the Hartree correction 



to the grand canonical potential is bounded by Ail < 
SN 2 /2g(e F ). Evaluating the elementary sums over peri- 
odic orbits, we find that the average interaction correc- 
tion (AO) is small compared to the mesoscopic oscilla- 
tions of CI: 



(An) 

rms Sfl 



< 



1.36791 

k F D* ' 



(14) 



where k F D* > 4.81 for a contact with nonzero conduc- 
tance. This result justifies the use of the independent- 
electron approximation [4-6,8-10]. 

In conclusion, we have shown that trace formulae a la 
Gutzwiller converge and give quantitatively accurate re- 
sults for the equilibrium quantum fluctuations in point 
contacts and nanowires. Using this approach, we have 
shown that the cohesive force of a metallic nanocon- 
tact, modeled as a hard- wall constriction in an electron 
gas, exhibits universal mesoscopic oscillations whose size 
rmsSF ~ e F /\ F is independent of the conductance and 
shape of the contact, and depends only on a dimension- 
less parameter a characterizing the degree of quantum 
chaos. Our prediction of universality is consistent with 
all experiments performed to date [3] . 
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